A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

02.05.1; 02.09.1; 02.14.1 06.1.2; 06.09.1 



ASTRONOMY 

AND 
ASTROPHYSICS 



New type astrophysical solution to the solar neutrino problems and 
its predictions to the SNO 



o 
o 
o 

(N ■ 

a ■ 



Attila Grandpierre 

Konkoly Observatory of the Hungarian Academy of Sciences, H-1525 Budapest, Hungary; grandp@konkoly.hu 



(N 

vn 
o 

o . 
o ■ 

Oh, 

O ■ 

■ 



Received . . . / Accepted . 



Abstract. The anomalously slow rotation of the solar core is 
just one from a remarkable lists of fundamental indications 
showing that the solar core is somehow coupled to the sur- 
face activity cycle. On the other hand, the atmospheric, LSND 
and solar neutrino problems are not consistent, therefore one or 
more of the neutrino experiments must be attributed - at least 
in part - to some phenomena other than neutrino oscillations, 
or a fourth neutrino is required. We use the luminosity con- 
straint in a general case not considered yet, taking into account 
nuclear reactions occurring in small hot bubbles indicated to 
be present in the solar core (Grandpierre, 2000). The dynamic 
solar model fit the observed neutrino fluxes within one-and- 
half sigma even without oscillations. This fit is shown to be 
better than the present day MSW and VAC fits. An attracting 
perspective is obtained for a simultaneous solution of the neu- 
trino problems and the solar core-related astrophysical prob- 
lems. Predictions of the dynamic solar model are presented for 
the SNO measurements. 
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1. Introduction 

Two different approaches seem to be promising in exploring 
the roots of the solar neutrino problems. One is the traditional 
particle physics approach, suggesting that neutrino oscillations 
are responsible for the missing solar neutrinos. Unfortunately, 
the small mixing angle, large mixing angle and vacuum oscilla- 
tion solutions all suffer from being poor fits to the observations. 
Eleven years ago Bahcall (1989) wrote in his book, that the 
MSW neutrino oscillation solution of the solar neutrino prob- 
lem (SNP) is attractive since the mixing angles and mass dif- 
ferences can each vary by orders of magnitude, and so it does 
not need fine tuning as the vacuum solution does. But not so 
much later Paterno and Scalia (1994) noticed that the allowed 
region shrinkened to a point-like area, and so "hypotheses on 
non-conventional neutrino properties are strongly disfavoured, 
except for the matter neutrino oscillations, the latter surviving 
within very narrow limits". Now even this remained point-like 
area (and the surrounding small 95 % C.L. region, which ex- 



tends only due to the theoretical and experimental errors) do 
not have a high allowance. The probability of these solutions 
decreased from 100 % to cca. 10% for the widely regarded best 
fit of small mixing angle (SMA) solution (and similarly for the 
LMA and VAC solutions), when the SuperKamiokande (SK) 
rates, spectral and day-night effect data are also taken into ac- 
count. 

One can evaluate the acceptability of the presented so- 
lutions by the probability P(x^) belonging to the calculated 
Xmin- The rule of thumb telling that "the value of for a 
"moderately" good fit is ~ d.o.f." (Press et al., 1992) may 
be useful. Bahcall, Krastev and Smirnov (1998) obtained for 
the SMA global fit Xmin = 26.5/l7d.o.f. which is acceptable 
at the probabiHty P{x'^) = 7 % C.L. Suzuki (1998) obtained 
Xmin — 50.2/31(i.o./. for the case when oscillations are not 
allowed. The corresponding probability is P(x^) = 1.6%. The 
best fit he obtained is for the vacuum oscillation (VAC), with 
X^.^ = 38.7/31 d.o.f., with P{x^) « 20%. This latter case 
we recognize also as a solution worse than moderately good, 
since Xmin > d.o.f. At the same time, Bahcall, Krastev and 
Smirnov (1998) obtained a global fit (when the constraints from 
the rates, the spectrum shape and the Day-Night asymmetry 
are all included) to for the best fit vacuum oscillation solution 
Xmin = 28.4 for 18 d.o.f., which is acceptable only at 6% C.L. 

It may be disturbing that two different and contrasting con- 
fidences are introduced to characterise the acceptance of a fit. 
One is P{x'^), the probability of a fit. The other one, P(Ax^), 
characterises how far the (allegedly normal) distribution of the 
experimental and theoretical errors extends. Since the larger 
errors leads to larger allowed regions, and the farther iso- 
probability contours belong to increasing n characterising the 
distance as na, therefore the higher is this second probability, 
the worse is the fit. Therefore, the higher is the former prob- 
ability -P(x^), the better is the fit; while, on the contrary, the 
higher is the latter probability, P(Ax^), the worse is the fit. 

In Gonzalez-Garcia and Pena-Garay (2000) Table 3, the 
larger C.L. belongs to a higher x^ /d.o.f and so it shows a 
poorer fit. They found for the SMA a C.L. 83%, which tells 
us that it is somewhere outside the 1 .5 fi region around the best 
fit which has a low (cca. 15 %) probability itself. Langacker 
(1999) noted that the LMA solution and the no oscillation hy- 
pothesis has the same rate of goodness, i.e. both of them is 
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disfavored at the 95-99 % CL level, while the LMA solution is 
also a very poor fit, although it is allowed at 95% CL. In other 
terms, this means that SMA is refuted with more than 1.5 cr. 

Recently Maris and Petcov (2000) found that "the conser- 
vative SMA around the point Am^ = 4 x 10~^, sm^29 = 
0.0085 is ruled out at 1.5 sigma" by the degree of the day-night 
effect observed by the SK. Now Haxton (2000) noted that "One 
puzzling aspect of atmospheric, solar, and LSND neutrino re- 
sults is that they require three independent Sm^s. That is, they 
do not respect the relation 5m\i + Sm^i + 5m}^2 = 0' thus ei- 
ther one or more of the neutrino experiments must be attributed 
to some phenomena other than neutrino oscillations, or a fourth 
neutrino is required." 

We learned that the neutrino oscillation solutions of the 
SNP actually are poor fits. Moreover, we recognise that even 
if they would be moderately good fits, it would be still an ur- 
gent need to find new physics outside from the neutrino oscilla- 
tions to explain the results of the neutrino experiments. In this 
Letter we suggest to look for the new physics in a conserva- 
tive region, in the field of astrophysics. The reason is the (not 
yet recognised) significance of the facts corresponding to the 
coupling of the solar core to the surface activity phenomena 
(one can find the presentation of the astrophysical core-related 
problems in Grandpierre, 1996, 1999, 2000). It is known, that 
the solar core rotates so slowly that its rotation rate is the sur- 
face rate ±30% (Elsworth et al., 1995). At the same time, the 
so-called 'best solar model' of Pinsonneault et al. (1989) pre- 
dicted a rotation rate of the core which is 4-15 times the surface 
rate. The discrepancy may be surmounted only when allowing a 
coupling between the solar core and the surface active regions. 
But if the core may participate in the surface activity cycle, the 
energy production has to be touched, and so the neutrino pro- 
duction may be influenced as well. In this way we found a way 
how the new physics arises in the field of astrophysics. 

The objections raised against a possible astrophysical solu- 
tion to the solar neutrino problems (see e.g. Bludman, Hata and 
Langacker 1994) are valid only for non-standard models with- 
out new physics (without changing the spectral shape of the 
individual neutrino fluxes). But it is known that the solution of 
the SNP needs new physics (Bludman, Hata and Langacker, 
1994). The often used term "model-independence" refers to 
general models in which the individual pp, Be^, CNO and 
neutrino fluxes are allowed to vary as free parameters. But if 
new physics is present, another kind of nuclear reactions may 
also contribute to the energy and neutrino production of the 
Sun. In this paper we attempt to show how the presence of 
a high- temperature energy source (Grandpierre, 1996, 2000) 
could contribute to the neutrino detector data. 

It is a general view that the principal neutrino sources are 
pp, pep, '^Be, ^B, ^^N, and ^^O. As BahcaU and Krastev 
(1996) remarked, this picture has become so accepted that it 
is sometimes referred to as "model-independent". Neverthe- 
less, if other neutrino sources do exist, then these previously 
thought "model-independent" models may all fail at a cer- 
tain rate. Therefore, the related luminosity constraint also may 
prove insufficient and over-constrained when compared to the 



actual Sun, since the hot bubbles may produce a smaller or 
larger part of the solar luminosity. The luminosity constraint 
expresses the fact that the energy productions is related to the 
neutrino production since both is produced by the nuclear re- 
actions: 

L0/47ri?2 = ^(Q/2- < E >„)</.(«), (1) 

a. 

where R=l A. U. (1.469x10^^ cm), < >„ is flie aver- 
age neutrino energy and <^(q;) is the ath flux (Q;=pp, pep, 
^-Be,^ B,...), and Q is the energy released in the fusion reaction 
4p -|- 2e — > a -|- 2u. The most general luminosity constraint we 
found in the literature is presented in Minakata and Nunokawa, 
1999 in the following form: 

Lq/A-kR^ = 13.10(pp) + 11.920(pep) + 12.5(/)(^Be) + 
6.660(^5) + 3.46(/)(^37V) + 2l.b7(l){^^0) 

+2M(p{^'^F) + lQ.n(j){hep). (2) 

We found it more suitable to convert these neutrino fluxes to 
fractional fluxes normalised to the BP98 SSM fluxes: 

1 = 0.9119$(pp) + 0.001966$(j5ep) + 0.0758$(Be) 

-F5.193 X 10-5$(B) -I- 3 X lQ-^^{CNO). (3) 

Neglecting the small terms, and including the pep into the pp, 
the CNO into the Be term, 

0.99 = 0.914$i -I- 0.076$7, (4) 

where $i refers to the pp-npep fluxes, and $7 to the Be-nCNO 
fluxes (where we modified the left side with the condition that 
for $i=$7=l the equation must be valid). 

We used the fractional neutrino fluxes following Minakata 
and Nunokawa (1998). The chlorine equation is: 

2.56= 1.8$7 + 5.9$8- (5) 
The gallium-equation is: 

72.4 = 69.6$i + 46.9$7 + 12.4$8- (6) 

It is easy to derive from these equations the lines of the 
Fig. 1 of Fiorentini and Ricci (1998), shown here as Fig. 1. It is 
confirmed that the non standard cool Sun models do not give an 
acceptable fit, as they are farther than 3cr from the intersection 
of the acceptance zones. 

Now we are prepared to consider the case when hot bub- 
bles are present. If the energy production of the hot bubbles is 
not negligible when compared to the total solar luminosity, we 
should involve it into the luminosity constraint. Regarding the 
energy budget of the hot bubbles, it depends on their temper- 
ature, which is indicated to be in the range of 10^ to 10^^ K 
(Grandpierre, 1996, 2000). One can think that the bubbles may 
produce energy through the hot CNO cycle, triple alpha cy- 
cle (this cycle does not produce neutrinos) and other nova-type 
nuclear reactions (Audouze, Truran, Zimmerman, 1973). From 
these reactions the one which may produce the largest number 
of events in the neutrino detectors may be the hot CNO cycle. 
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Fig. 1. The and ^Be+CNO neutrino fluxes, consistent with 
the luminosity constraint and experimental results for standard 
neutrinos. The SSM flux region is around the (1,1) point. The 
dashed (sohd) lines correspond to the central (ilex) experi- 
mental values for CI, Ga and u — e scattering experiments. 
The crosses indicate the behaviour of non standard solar mod- 
els with low central temperature. 



Fig. 2. The same as in Fig. 1, except that the luminosity con- 
straint is now written for the quiet solar core of the dynamical 
solar model. The hot bubbles are allowed to produce mu and 
tau neutrinos, therefore the Kam line shifts leftwards. The tri- 
angle LQ'99 refers to the calculations of Lavagno and Quarati, 
1999. The figure shows that LQ'99 does not fit exactly to the 
low-temperature parabola. 



Writing the luminosity constraint for the hot bubble separately, 
for the case when the CNO cycle gives ten percent of the total 
solar luminosity as an upper limit, from (1) we obtain: 

O.IL© = 12.b2b^''{CNO). (7) 

From this constraint, we can derive an upper limit for the bub- 
ble CNO neutrino flux, 

^^{CNO) < 6 X lO^cm-^s-^, (8) 

which is close to the value of the SSM beryllium-neutrino 
flux. This value is compatible with the present-day global con- 
straints 0.0 < <i>CBe) < 6.35 (Bahcall and Krastev, 1996). 

Now we can calculate how Fig. 1 is modified when the hot 
bubbles are producing 22% of the total solar luminosity. In this 
case the luminosity constraint for the quiet solar core without 
the bubbles should be formulated as: 

0.77 = 0.914$i + 0.076$7- (9) 

Using this constraint in the Ga-equation (6), we observe that 
the effect of the hot bubble energy generation, or, more pre- 
cisely, the constraint that the quiet solar core should produce 
a less than total solar luminosity, is to shift the Ga-zone up- 
wards. This effect is helpful in obtaining better and physical 



(i.e. ^{Be) > 0) fluxes, and so to resolve one of the solar neu- 
trino problems, the problem of beryllium flux. 

Now one may consider what happens with the Kamiokande 
zone if hot bubbles are present. Since the bubbles may reach 
very high temperatures, they may produce mu and tau neutri- 
nos, which can be observed by the Kamiokande and not by the 
other detectors. Therefore, since Fig. 2 describes the quiet so- 
lar core only, the Kam-zone may be shifted to the left. Actually, 
the rate of the necessary shift is determined by the intersection 
of the CI and Ga zones. We observe from Fig. 2 that the shift 
occurs around $8=0.26, which means that the bubbles has to 
produce a j/^ neutral current with a contribution to the SK 

= 0.21. (10) 

We note here, that since the hot bubbles may contribute also 
to the high-energy excess observed at the SK, therefore the 
amount of ^ should not be so large. The Cl-zone of Fig. 
1 will not be modified when going to Fig. 2. 

2. Discussion and Conclusions 

It is interesting that the dynamic solar model (DSM) obtained 

with the inclusion of the hot bubbles into the standard so- 
lar model (SSM) modifies the Ga, CI, Kam and non-standard 
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zones in a way the create shifting an overlap region in a physi- 
cal range ($7 > 0. Remarkably, the cool Sun model also over- 
laps with this overlap region. For a bubble-luminosity around 
0.22 (of the total solar luminosity), they all fit with cca. 1.2 a. 
Fig. 2 shows that the three kinds of neutrino detectors actually 
do not contradict to the standard neutrino picture. The situation 
is that the combined neutrino results actually indicate a lower 
than standard central temperature of the Sun with a remarkable 
confidence. Fig. 2 gives a better fit (cca. 1.2 a) than that of the 
presently favoured SMA solution (higher than 1.5ct). The best 
fit may be reached when both the bubbles and the oscillations 
are taken into account. 

One can observe from Fig. 2 that the new astrophysics may 
offer powerful perspectives to influence the solution of the SNP. 
It seems to be not true, that we know the basic physics of the 
Sun enough. Even if one can consider the presence of the hot 
bubbles in the solar core as yet not established, the frequent 
statements that astrophysical solutions are ruled out, prove to 
be unfounded in the here presented more general basis. 

In Fig. 2 the non standard low- temperature solutions are 
also presented. Although in Fig. 1 they are found far from the 
1 a zones, in the case of Fig. 2 they are around 1.2 cr from 
the overlap region when = 0.942. This value of central 
temperature is quite consistent with the 0.78 solar luminosity 
produced by the quiet solar core. 

We do not attempt to suggest that the hot bubbles perma- 
nently produce a significant part of the solar luminosity. Our 
attempt is more narrowly confined: it is to point out that the 
bubbles in principle may create a situation in which the astro- 
physical solution may be alone enough to solve the problem of 
the missing solar neutrinos. It is quite plausible, that the elec- 
tron neutrinos do oscillate, and this phenomenon is responsi- 
ble for a large part of the SNP and other apparent "neutrino 
anomalies". Nevertheless, the point is that if we ignore the role 
played by the dynamics of the solar core in the production of 
solar neutrinos, we may found ourselves in the uncomfortable 
situation of poor fits of the MS parameters and inconsistent 
consequences of atmospheric, LSND and solar neutrino prob- 
lems. 

To resolve the astrophysical problems of the solar core it is 
not necessary that the hot bubbles produce a significant part of 
the solar luminosity. Actually, they may produce only a neg- 
ligible part of the solar luminosity and they could prove still 
be able to trigger and influence the surface solar activity. Gor- 
batzky (1964) already calculated that hot bubbles arising from 
point explosions with an initial surplus energy around 10^'' ergs 
may be able to reach the stellar surfaces from 0. 1 stellar radius. 
He ignored completely the energy production of these hot bub- 
bles with temperatures T > 10^ K. But even a slight influence 
of the solar core dynamics might be able to improve signifi- 
cantly the fit and consistency of the neutrino problems. 

The result red from Fig. 2 has predictive value for the fu- 
ture neutrino detectors. It is easy to estimate the consequences 
of this dynamic solar model (for standard neutrinos) regarding 
the SNO observations. Our picture modifies the conclusion that 
the [NC]/[CC] is larger than unity ([i]=observed rate/standard 



solar model rate, see Bahcall, Krastev, Smirnov, 2000) will def- 
initely indicate the presence of neutrino observations. Since 
at high temperatures like 10^° — lO^^i^ the hot bubbles may 
produce mu and tau neutrinos independently of the presence 
of neutrino oscillations, they may contribute to [NC] increas- 
ing it. For example, if we take T ~ 0.942, than we will have 
Rfj^ = 0.23, and so [NC]/[CC]=2.04. If we aflow a value of 
T closer to the SSM value 1, the ratio {[NC]/[CC])dsm de- 
creases towards unity. This conjecture of the DSM should be 
taken into account at the interpretations of the SNO and SK 
observations. 

Regarding the helioseismic context, we note that sound 
speed is not sensitive to the nuclear reactions. Bahcall and 
Ulrich (1988) and Basu, Pinsonneault and Bahcall (2000) re- 
marked that even when switching out the He^ + He^ reaction 
(producing 14 % of the total solar luminosity), the sound speed 
differ only by 0. 1 % from the sound speed obtained from the 
standard solar model. On the other hand, the calculations of the 
dynamic solar model show that when bubbles and neutrino os- 
cillations are both present the SMA and LMA regions (of the 
purely MSW solutions) shift significantly even for such a small 
change as T 0.995. Moreover, the many-body effects of the 
particles in the nuclear reactions led to solar models which may 
compensate the cool solar models towards an SSM central tem- 
perature (Lavagno, Quarati 1999). 

References 

Audouze, J., Truran,;., & Zimmerman, B.A. 1973, ApJ, 184, 493 
Bahcall, J. N., Krastev, P. I., Smirnov, A. Yu. 1998, Phys. Rev. D 58, 
096016 

Bahcall, J. N., Krastev. P. I., Smirnov, A. Yu. 2000, a rXiv:hep- 
ph/0002293 

Bahcall, J. N. 1989, m Neutrino Astrophysics, Cambridge Univ. Press, 
246 

Bahcall, J. N. and Krsatev, P I. 1996, Phys. Rev. D 53, 421 1 
Bahcall, J. N. and Ulrich, R. K. 1988, Rev. Mod. Phys. 60, 297 
Basu, S., Pinsonneault, M. H., Bahcall, J. N. 2000, ApJ 529, 1084 
Bludman, S., Hata, N. and Langacker, P 1994, Phys. Rev. D. 49, 3622 
Elsworth, Y., Howe, R., Isaak, G. R., McLeod, C. P Miller, B. A., 

Wheeler, S. J.,New, R., and Gough, D. O. 1995, Nat 376, 669 
Fiorentini, G. And Ricci, B. 1998, in Proc. of Intern. Workshop "Be- 
yond t he standard model : from theory to experiment ", Valencia, 



Spain, 3stro-ph/9801185 



Gonzalez-Garcia, M. C. and Pena-Garay, C. 2000, |hep-ph/0001015 

Gorbatzky, V. G. 1964, AZh 41, 53 

Grandpierre, A. 1996, A & A 308, 199 

Grandpierre, A. 1999, A & A 348, 993 

Grandpierre, A. 200 0, Nucl. Phys. B (Proc. S uppl.), in print 

Haxton, W. C. 2000, kXiv:nucl-th/0004052| 



Langacker, P 1999, Nucl. Phys. P roc. Suppl. 77. 241 
Lavagno, A. and Quarati. P 1999, |arXiv:nucl-th/9912064 



Maris, M. and Petcov. S. T. 2000. arXiv:hep-ph/0004151 



Minakata. H. and Nunokawa, H. 1999, Phys. Rev. D 59. 073004 
Patemo. L. and Scalia, A. 1994, in Advances in Solar Physics, eds. 

Belvedere, G., Rodono, M. and Simnett, G. M., 41 
Pinsonneault, M. H., Kawaler, S. D., Sofia, SA. and Demarque, P. 

1989, ApJ 406,714 



Attila Grandpierre: New type astrophysical solution to the solar neutrino problems and its predictions to the SNO 



5 



Press, W.H., Teukolsky, S. A., Vetterling, W. T. and Flannery, B. P. 

1992, in Numerical Recipes in FORTRAN, 2nd ed., Cambridge 

Univ. Oress, 655 
Suzuki, Y. 1999, Nucl. Phys. B (Proc. Suppl.) 77, 35 



